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This article investigates the signature of the seventeen multi-connected flat spaces in cosmic 
^ \ microwave background (CMB) maps. For each such space it recalls a fundamental domain and a set 

"nI . of generating matrices, and then goes on to find an orthonormal basis for the set of eigenmodes of 

the Laplace operator on that space. The basis eigenmodes are expressed as linear combinations of 
eigenmodes of the simply connected Euclidean space. A preceding work, which provides a general 
method for implementing multi-connected topologies in standard CMB codes, is then applied to 
simulate CMB maps and angular power spectra for each space. Unlike in the 3-torus, the results 
in most multi-connected flat spaces depend on the location of the observer. This effect is discussed 
in detail. In particular, it is shown that the correlated circles on a CMB map are generically not 
back- to-back, so that negative search of back-to-back circles in the WMAP data does not exclude a 
^1^' vast majority of flat or nearly flat topologies. 
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PACS numbers: 98.80.-k, 04.20.-q, 02.40.Pc 



I. INTRODUCTION 



o 

> 

' Among all multi-connected three-dimensional spaces, "flat spaces"^ have been studied the most extensively in the 
cosmological context. This is due to the computational simplicity of the simplest compact flat three-manifold, the 
3-torus, which has been used extensively in numerical simulations. The main goal of this article is to provide tools to 
compute the CMB properties and produce high resolution CMB maps for all seventeen multi-connected flat spaces ^, 
following the general method introduced in our preceding work Q . 

Recent measurements show that the density parameter VLq is close to unity and the observable universe is approx- 
imately flat. CMB data obtained by the Archeops balloon experiments 2] and more recently by WMAP 3] place 
strong constraints on the curvature. In addition, WMAP 4] and later the Planck satellite 5] do and will provide 
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^ In this article, we follow the cosmological use and we call "flat spaces" 
the simply connected universal cover E^. 
Test maps for these spaces are available on demand. 



the eighteen types of Euclidean spaces, and "Euclidean space" 
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full sky maps of CMB anisotropics, offering an opportunity to probe the topological properties of our universe. This 
observational constraint on the curvature radius of the universe motivates the detailed study of flat spaces even though 
spherical spaces are also promising candidates 0, S S S ^ • 

At present, the status of the constraints on the topology of fiat spaces is evolving rapidly driven by the release of 
the WMAP data. Previous analysis, based on the COBE data, mainly constrained the topology of a 3-torus (see 
Refs. [illilllllillillillillilia Refs. jHIHiBl for reviews of different methods for searching for 

the topology). 

The WMAP data 3] possess some anomalies on large angular scales that may be explained by a topological 
structure. In particular, the quadrupole is abnormally low, the octopole is very planar and the alignment between the 
quadrupole and octopole is also anomalous . Besides many other potential explanations , it was suggested that 
a toroidal universe with a smaller dimension on the order of half the horizon scale may explain all these anomalies 26] 
but it was latter shown, on the basis of a finer statistical analysis, that it did not |23. Another topology was recently 
proposed to explain some of this anomaly in the case of a slightly positively curved space, namely the Poincare 
dodecahedral space [Tl| . 

The first results of the search for the topology through the "circles in the sky" method ^29] gave negative results 
for back-to-back or almost back-to-back circles 28, 30j. While the first applies only to back-to-back circle with no 
twist, the second study includes an arbitrary twist and conclude that "it rules out the possibility that we live in a 
universe with topology smaller than 24 Gpc" . As will be discussed in this paper, back-to-back circles are generic only 
for homogeneous topologies such as e.g. 3-tori and a subclass of lens spaces. In non-homogeneous spaces the relative 
position of the circles depends on the position of the observer in the fundamental polyhedron. 

In conclusion, as demonstrated by these preliminary results, only the toroidal spaces have been really constrained 
|3Q( . Besides, a series of studies have pointed out a departure of the WMAP data from statistical isotropy. Copi et 
al. |31| recently argued in particular that they are inconsistent with an isotropic Gaussian distribution at 98.8% 
confidence level. Previous studies pointed toward a possible North-South asymmetry of the data j^^^lSl- Spaces 
with non-trivial topology are a class of models in which global isotropy (and possibly global homogeneity) is broken. 
Simulated CMB maps of these spaces may help to construct estimators for quantifying the departure of the temperature 
distribution from isotropy, and also give a deeper understanding of recent results. 

Let us emphasize that in the case where the topological scale is slightly larger than the size of the observable universe, 
no matching circles will be observed. This might also happen for a configuration where the circles would all lie in the 
direction of the galactic disk where the signal-to-noise ratio might be too low. Contrary to the simply connected case, 
the correlation matrix, Cf^ = {a£mO^£'m')^ of the coefficients of the development of the temperature fluctuations 
on spherical harmonics, will not be proportional to SwSmm'' The study of this correlation matrix could offer the 
possibility to probe topology (slightly) beyond the horizon. Computing the correlation matrix Cj^ for different 
multi-connected spaces will help design the best strategy to constrain the deviation from the simply connected case, 
and gives a concrete example of cosmological models in which the global homogeneity and isotropy are broken. 

As described in detail in our preceding work Jj, what is needed for any CMB computation are the eigenmodes of 
the Laplacian 

ATfU-fc^Tfl, (1) 

with boundary conditions compatible with the given topology. These eigenmodes can be developed on the basis 
of the (spherically symmetric) eigenmodes of the universal covering space as 

oo i 

^fcs"^^ ^ ^klmykim, (2) 

fri s 

so that all the topological information is encoded in the coefficients where s labels the various eigenmodes 

sharing the same eigenvalue — /c^. Ref. ^ computes these coefficients for the torus and lens spaces and Refs. fA l3^ 
discuss more general cases. 

To summarize, this article aims at several goals. First, it will give the complete classification of flat spaces and 
the exact form of the eigenmodes of the Laplacian for each of them. It will also provide a set of simulated CMB 
maps for most of these spaces. Among other effects, it will illustrate the effect of non-compact directions and discuss 
the influence of the position of the observer in the case of non-homogeneous spaces, which has never been discussed 
before. It also explains the structure of the observed CMB spectrum in the case of a very anisotropic (i.e., flattened 
or elongated in one direction) fundamental domain. 

This article is organized as follows. We start by recalling the properties of the eighteen flat spaces (Section^ as 
well as the eigenmodes of the simply connected three-dimensional Euclidean space (Section llll|) . and in particular 
how to convert planar waves, which suit the description of topology, to spherical waves, which are more convenient 
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TABLE I: Classification of the 18 three-dimensional flat spaces. 



for CMB computation. Then, in Section Hvi we explain how to extract the modes of a given multi-connected space 
from the modes of E^. This method is then applied to give the eigenmodes of the ten compact flat spaces (Section IV)) . 
the five multi-connected flat spaces with two compact directions ( "chimney spaces" , Section IVI|) and the two multi- 
connected flat spaces with only one compact direction ("slab spaces", Section IVlT)) . Applying the general formalism 
developed in our previous work Q, we produce CMB maps for some of these spaces. With three exceptions the 
manifolds are not homogeneous, in the sense that a given manifold does not look the same from all points. To discuss 
the implication of the observed CMB and the genericity of the maps, we detail in Section IIXI the influence of the 
position of the observer on the form of the eigenmodes and we study its consequences on the observed CMB maps. 
We show in particular that the matched circles are generically not back-to-back, but their relative position depends 
on the topology, the precise shape of the fundamental domain, and the position of the observer. 

Notation 

The local geometry of the universe is described by a locally Euclidean Friedmann-Lemaitre metric 

= -c'dt' + it) [dx^ + x'dio^] , (3) 
where a{t) is the scale factor, t the cosmic time, do;^ = dO'^ + sin^ dcp'^ the infinitesimal solid angle. 

II. THE EIGHTEEN FLAT SPACES 

Let us start by recalling the list of flat spaces. They are obtained as the quotient E^/F of three-dimensional 
Euclidean space E^ by a g:roup F of symmetries of E^ that is discrete and fixed point free. The classification of such 
spaces has long been known 35, ^3^, motivated by the study of crystallography and completed in 1934 ^3?\. The 
ten compact flat spaces are quotients of the 3-torus; six are orientable and four are non-orient able. Fig. shows 
fundamental polyhedra. The non-compact spaces form two families, the chimney space and its quotients having 
two compact directions (Fig. ^ and the slab space and its quotient having only one compact direction (Fig. 
The terms slab space and chimney space were coined by Adams and Shapiro in their beautiful exposition of the flat 
three-dimensional topologies j3^]. Table ffl summarizes the properties of the whole family of flat spaces. 
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Klein Space Klein Space 

with Horizontal Flip 




Klein Space Klein Space 

with Vertical Flip with Half Turn 



FIG. 1: Fundamental domains for the compact flat three-manifolds. The unmarked walls are glued straight across. Courtesy 
of Adam Weeks Marano (first published in Ref. ^^). 

III. EIGENMODES OF 

The eigenmodes of Euclidean space admit two different bases: a basis of planar waves and a basis of spherical 
waves. The former are more convenient when seeking eigenbases for multi-connected spaces, while the latter are more 
convenient for simulating CMB maps. This section considers both bases and the conversion between them, as also 
detailed in the particular case of the torus in Ref. JLj . 



A. Planar waves 

Each vector k defines a planar wave 

Tk(x) = e*'^-'^. (4) 

The defining vector k, called the wave vector^ lives in the dual space, so the dot product k • x is always dimensionless. 
These modes are indeed not square integrable and are normalized as 

Tk(x)Tj;,(x)^ = 5^{k - k'). (5) 



Chimney Space Chimney Space with 

Half Turn 




Chimney Space with 
Vertical Flip 





Chimney Space with 
Horizontal Flip 



Chimney Space with 
Half Turn and Flip 



FIG. 2: The chimney space is made from an infinitely tah rectangular chimney with front and back faces (resp. left and right 
faces) glued straight across. The four variations on the chimney space space glue the front face to the back as indicated by the 
doors. In all variations except the last the left and right faces are glued straight across; in the last variation they are glued 
with a top- to-bottom flip so that the windows match. (Courtesy of Adam Weeks Marano) 





Slab Space 



Slab Space with Flip 



FIG. 3: The slab space is made from an infinitely tall and wide slab of space with its front face glued to its back face straight 
across. The variation glues the faces with a flip. (Courtesy of Adam Weeks Marano) 



B. Spherical waves 

Each spherical wave factors into a radial part and an angular part, 



ykim{x, ^. ^) = y f (27r)'/' 3t{kx) Ytie, (6) 



where (x, ^, ^) are the usual spherical coordinates 

X = X sin^ cosip 
y = X sin^ sincp 

z = X cos 6. (7) 

The radial factor ji{kx) is the spherical Bessel function of index ^, and the angular factor YJ^{0^(p) is the standard 
spherical harmonic. The mode yk£m is not square integrable and is normalized according to 

/ ^) x^dxdcos6>d(/9 1 ^ . /nx 

J^^ yMmyk'l'm' ^ 1^ [k-k )dn'Omm', (8) 

which is analogous to the normalization Q and which determines the numerical coefficient \/2/ti. 
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C. Conversion 

Subsequent sections will find explicit bases for the eigenmodes of multi-connected flat three-manifolds as linear 
combinations of planar waves. The planar waves may easily be converted to spherical waves using Eqns 5.17.3.14 and 
5.17.2.9 of Ref. 

Tk(x) = e^^-^ 

= ^ jfXk |x|) {2e + 1) PKcos0k,x) 
e=o 



oo i 

J2 E ^'^r*(k) [47ri,(fc|x|)rr(x)] 

i=0 m=-i 
oo £ 

J2 E {i'Yr{^))ykiM, 



(9) 



=0 m=-£ 



where k = |k|, k = k/|k|, and x = x/|x|. 

In particular, the conversion formula © lets one easily translate a planar wave Tk to the framework we developed 
in Ref. 0, which expresses each basis eigenmode as a sum of spherical waves 



(10) 



£=0 m=-£ 



where 5 indexes the different Tk whose wave vectors k share the same modulus k. In the Euclidean case the index 
may be chosen to be 5 = k. Comparison with Q immediately gives the coefficients 



(11) 



IV. EIGENMODES OF MULTI-CONNECTED SPACES 

A multi-connected fiat three- manifold is the quotient E^/F of Euclidean space under the action of a group T 
of isometrics. The group T is called the holonomy group and is always discrete and fixed point free. Each eigenmode 
T of the multi-connected space E^/F lifts to a F-periodic eigenmode T of E^, that is, to an eigenmode of E^ that is 
invariant under the action of the holonomy group F. Common practice blurs the distinction between eigenmodes of 
E^/F and F-periodic eigenmodes of E^, and we follow that practice here. Thus the task of finding the eigenmodes 
of the multi-connected space E^/F becomes the task of finding the F-periodic eigenmodes of E^. In this section we 
investigate how an isometry 7 G F acts on the space of eigenmodes. The two lemmas we obtain will make it easy to 
determine the eigenmodes of specific multi-connected spaces in subsequent sections. 

Every isometry 7 of Euclidean space E^ factors as a rotation/reflection followed by a translation. If we write the 
rotation/reflection as a 3 x 3 matrix M in the orthogonal group 0(3) and write the translation as a vector T, then 7 
acts on E^ as 



(12) 
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This isometry of E^ induces a natural action Tk(x) Tk(Mx - 



Lemma 1 (Action Lemma). The natural action of an isometry 7 0/ E^ 
planar eigenmode e*^"'^TkM- 



T) on the space of eigenmodes. 

takes a planar eigenmode Tk to another 
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Proof. Keeping in mind that k is a row vector while x is a column vector, the proof is an easy computation: 



Tk(x) 



^ gik-(Mx+T) 



= e^^-^TkM(x). (13) 

Q.E.D. 

Lemma 2 (Invariance Lemma). If j is an isometry o/E^ with matrix part M and translational part T , the mode 
Tk is a planar wave, and n is the smallest positive integer such that k = kM"' (typically n is simply the order of the 
matrix M ), then the action of 7 

1. preserves the n-dimensional space of eigenmodes spanned by {Tk, TkM, • • • , TkM^-i} ^ ^^t, and 

2. fixes a specific element 

aoTk + aiTkM + • • • + ^n-iTkM— 1, (14) 

if and only if for each j (mod n) 

=e^^^''^a,. (15) 
Proof. Both parts are immediate corollaries of Lemma 1. Specifically, the action of 7 takes a linear combination 

<^oTk + <^iTkM + . . . + an-2TkM'^-2 + a^-iTkAf^-i^ (16) 

to 

ttoe IkM+^^ie lkM2 + . . . + fin-2e IkM— i+«n-ie Ik, (I'j 

SO it's clear that the n-dimensional subspace spanned by {Tk, TkM, • • • , TkM^-i} is preserved as a set. Equating 
(fTB|) to ([T7)) and comparing coefficients proves the second part. Q.E.D. 



V. COMPACT FLAT THREE-MANIFOLDS 



We will ffist find the eigenmodes of the 3-torus, and then use them to find the eigenmodes of the remaining compact 
flat three- manifolds. 



A. 3-Torus 



The 3-torus is the quotient of Euclidean space E^ under the action of three linearly independent translations Ti, 
T2 and T3. Its fundamental domain is a parallelepiped. Its eigenmodes are the eigenmodes of E^ invariant under the 
translations Ti, T2 and T3 (recall from Section |^ the convention that eigenmodes of the quotient are represented 
as periodic eigenmodes of E^). The Invariance Lemma (with n = 1) implies that an eigenmode Tk of E^ is invariant 
under the translation Ti if and only if e**^ '^^ = 1 which occurs precisely when k- Ti is an integer multiple of 27r. Thus 
geometrically the allowed values of the wave vector k form a family of parallel planes orthogonal to Ti. Similarly, 
the eigenmode Tk is invariant under the translation T2 (resp. T3) if and only if k lies on a family of parallel planes 
orthogonal to T2 (resp. T3), defined by k • T2 G 2tiZ (resp. k • T3 G 27rZ). The eigenmode Tk is invariant under all 
three translations Ti , T2 and T3 if and only if it lies on all three families of parallel planes simultaneously. The in- 
tersection of the three families forms a lattice of discrete points. Fig. ^illustrates the construction for the 2-torus; the 
construction for the 3-torus is analogous. This lattice of points defines the standard basis for the eigenspace of a torus. 

Definition. The standard basis for a 3-torus generated by three linearly independent translations Ti, T2 and 
T3 is the set B = {Tk | k • G 27rZ for i = 1, 2, 3}. 
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FIG. 4: An eigenmode Tk of is invariant under the translation if and only if k lies on a family of parallel lines orthogonal 
to Ti, defined by k • G 2tiZ. The mode Tk is invariant under both translations Ti and T2 if and only if it lies in the 
lattice of intersection points of the two families of parallel lines. The construction in three dimensions is similar, but with three 
families of planes instead of two families of lines. Strictly speaking we should not draw the parallel lines in the same space as 
the translation vectors because the wave vectors k live in the dual space (with units of length"^) while the live in the 
primary space (with units of length), but nevertheless it's visually helpful to do so. 




FIG. 5: In a rectangular 3-torus (left) the allowed wave vectors k form a rectangular lattice. In a hexagonal 3-torus (right) the 
lattice is hexagonal within each layer. 



The most important special case of a 3-torus is the rectangular 3-torus generated by three mutually orthogonal 
translations 

Ti = (L„0,0) 
T2 = (0,L^,0) 

T3 = (0,0, L,), (18) 

in which case the allowed wave vectors k take the form 

^^.fn n ny ^^^^ 

\ -LJx J^y J^z J 

for integer values of n^^, Uy and n^, thus forming a rectangular lattice (Fig. [51 left). 
The second most important special case is the hexagonal 3-torus generated by 




T4 = (0,0,L,). (20) 

These four generators and their inverses define a fundamental domain as a hexagonal prism. The first three generators 
are linearly dependent (Ti + T2 + T3 =0); eliminating any one of them suggests an alternative fundamental domain 
as a prism with a rhombic base. Even though the hexagonal and rhombic prisms look different, they define the same 
manifold. Either way, the allowed wave vectors k form a hexagonal lattice (Fig. [HI right) and may be parameterized 
as 



for integer values of ni, n2 and 713. 
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In the case of a general 3-torus, one writes the translations Ti, T2 and T3 as the columns of a 3 x 3 matrix T and 
solves the equation kT = 27r(ni, 77,2, ^3) to find the allowable wave vectors k = 27r(ni, 722, n3)T~^ for integers ni, 712 
and 713. 

When one wants to simulate CMB maps, one needs to know not only the modes themselves but also how the modes 
are paired under complex conjugation. The reason is that the cosmological fields are in fact real- valued stochastic 
variables. Any such field can be decomposed into Fourier modes as 

/d^k 
^^-^</,k(t)e^'^-ek, (22) 

where ek is a complex, usually Gaussian, random variable satisfying 

(ekej;,) = 5°(k-k'). (23) 

The evolution equations of the cosmological perturbations involve time derivatives and a Laplacian so that the 
coefficient ^k(^) can be decomposed as 

</'k(t) =<^fe(^)e'''^ (24) 

where is a phase that is constant throughout the evolution. By absorbing the phase into the random variable, we 
can always choose 0k to be a real function of k only, i.e., 0/c(t). Since e**^"^ and e~**^'^ are conjugate in the preceding 
decomposition, the fact that 0(x, t) is real implies that 

ek = e-k. (25) 

This latter constraint may not hold for all the other spaces studied in this article and we will need to give its analog 
for each case. 



B. Quotients of the 3-Torus 

For ease of illustration we first explain our general method for the two-dimensional Klein bottle. Figure El shows 
a portion of the Klein bottle's universal cover, in which alternate images of the fundamental domain appear mirror 
reversed. Half the holonomies are pure translations while the other half are glide refiections. In other words, the 
holonomy group F contains an index 2 subgroup F' C F comprising the pure translations. While E^/F gives the 
original Klein bottle, E^/F' gives a torus whose fundamental domain is the square formed by the solid lines in FigureEl 
(ignoring the dotted lines). According to the convention introduced in Section Hvi the Klein bottle's eigenmodes are 
represented as F-periodic functions on E^. But every F-periodic function is automatically a F'-periodic function as 
well, because F' is a subgroup of F. Thus every eigenmode of the Klein bottle is a priori an eigenmode of the 2-torus. 
The task in finding the eigenspace of the Klein bottle is to start with the eigenspace of the torus and find the subspace 
that is invariant under the glide refiection (the one taking the lower half of a square to the upper half). In practice 
this is quite simple. A rectangular torus has holonomy group F' generated by the two translations 

The standard eigenbasis for this torus takes the form (|19|) . namely B = {Tk | k = 27r(^, j^) for Ux^riy G Z}. To 
extend this F' to the full holonomy group F of the Klein bottle, we add the glide refiection 

(27) 

and ask which elements of the basis B it preserves. The Invar iance Lemma provides a ready answer: when kx the 
two-dimensional subspace {T^^^ky^^ -k^,ky} is preserved as a set (Part 1 of the Invariance Lemma) while the mode 
^k^,ky {—'^)^^^ -k^,ky is fixed exactly (Part 2 of the Invariance Lemma). In the exceptional case that kx = 0^ the 
one-dimensional subspace {To,/cy} is preserved as a set, while To,/cy is fixed if and only if Uy is even (because when 
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FIG. 6: The Klein bottle's holonomy group contains glide reflections as well as translations. The translations alone form an 
index 2 subgroup defining a torus. 



1, Part 2 of the Invariance Lemma requires ao = e ao which implies k • T/(27r) = nz/2 ^ Z). In summary, an 



orthonormal basis for the space of eigenmodes of the Klein bottle is the union of 



T 



for 
for 



Tlx e z^.Uy e z, 

Uy e 2Z. 



(28) 



Let us now apply this same method to each of the nine quotients of the 3-torus. 

1 . Half turn space 

The analysis of the half turn space closely follows that of the Klein bottle given immediately above. The only 
difference is that the Klein bottle's holonomy group contained translations and glide reflections, while the half turn 
space's holonomy group contains translations and corkscrew motions. Specifically, we begin with the generators for 
the holonomy group V of a rectangular 3-torus 



(x\ 

y 



y 
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(29) 



and add a generator for the half-turn corkscrew motion 
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(30) 



to get the full holonomy group T of the half turn space. The Invariance Lemma shows that when {kx^ ky) ^ (0, 0) the 
two-dimensional subspace {T/e^^/c^^/c^, T_/c^^_/c^^/c^} is preserved as a set (Part 1 of the Invariance Lemma) while the 
mode ^ka^,ky,kz + {~'^)^^^ -ka^,-ky,kz is fixed exactly (Part 2 of the Invariance Lemma). In the exceptional case that 
{kx^ ky) = (0, 0), the one-dimensional subspace {Tcc/c;^} is preserved as a set, while To,o,/c;s is fixed if and only if is 
even. In summary, an orthonormal basis for the space of eigenmodes of the half turn space, t|^^^[ ^ , is the union of 



1 

V2 



^27r(^ li. ) + {-lT'^27r(-:^ ^) ^^r (n^ e Z+, Uy, Uz e Z) or (n^ = 0, n^^ G Z+, G Z), 



-27r(0,0,J^) 



for Uz G 2Z. 



(31) 



In terms of the notations used in Ref. 1], it leads to the coefficients 

y^^*(k) + (-l)^-y^^*(kM)j for (ux e Z+, n^, Uz e Z) or {ux = 0, G Z+, G Z), 



^kim 



V2 

z^y^^*(k) 



for (na,,%) = (0,0), G 2Z. 



(32) 
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FIG. 7: The last scattering surface seen from outside for a half-turn space E2 with Lx = Ly — 0.64, Lz = 1.28 in units of the 
radius of the last scattering surface. Each row presents the two pairs of matching circles along the z direction, but seen from 
opposite directions. The left panels show circles which match with a tt rotation while the right panel shows circles which match 
without rotation. We recover the invariance by a translation of Lz due to the fact that the cubic torus of size Lz = 1.28 is the 
double cover of the half-turn space considered here. Only the Sachs- Wolfe contribution has been depicted here. 

k being given by Eq. (^HJ- Passing from the expression of the modes to the coefficients ^^^^ is straightforward and 
in the following we will give only the expressions of the modes. 

If desired, one could construct a more general half turn space from a right prism with a parallelogram base, instead 
of a rectangular box. 

To find the analog of Eq. ((23 , one simply needs to check that 

so that it follows that 

1. when A:^ 7^ 0, ek is a complex random variable satisfying 

eL,fe„fc, = {-^T'ek^,ky,-k., (34) 

2. when /c^ = 0, ek is a real random variable. 
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2. Quarter turn space 

The quarter turn space is similar to the half turn space, but with with a quarter turn corkscrew motion 







(o 


-1 






(x 
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(35) 



In particular this implies that Lx = Ly. The Invariance Lemma shows that when (kx^ky) ^ (O^O) the four- 
dimensional subspace {T/e^,/e^,fe, , Tky-k^,k,, T_fc, T-ky,k^,kA is preserved as a set, while the mode 
^k^,ky,k, + i''^Tky-k^,k, + {-1)''^T -k^-ky,k, + {-i)''^T-ky,k^,k, ^^^^ cxactly. In the exceptional case that 
(kx^ky) = (0,0), the one-dimensional subspace {To,o,/c^} is preserved as a set, while To,o,/c^ is fixed if and only if 
Uz is a multiple of 4. In summary, an orthonormal basis for the space of eigenmodes of the quarter turn space, 
'^l^^l -A- 5 is the union of 



1 

2 L 



"^^y Lx' Ly ^ Lz ^ ^^^Ly^ Lx'Lz^ "^^^ Lx' Ly^Lz^ 



for Ux e Z+, Uy e U {0}, riz G Z, 
T2^(o,o,5j) ^ 

As in the case of the half turn space, one can easily check that 

so that the analog of Eq. (|25|) is given by 

1. when /c^ 7^ 0, ek is a complex random variable satisfying 



(36) 



(37) 



(38) 



2. when /c^ = 0, ek is a real random variable. 



3. Third turn space 

The third turn space is a three-fold quotient of a hexagonal 3-torus, not a rectangular one. To the generators 
of the hexagonal 3-torus we add a one-third turn corkscrew motion 
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V 1 
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y + 
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(39) 



The eigenmodes Tk of the hexagonal 3-torus are already known from Eqn. ((2U (and illustrated in FigurelSl). Applying 
the Invariance Lemma to them with the additional generator ((S^ yields the eigenbasis, , 



^ [Tk + C^TkM + C'^^TkM^] for m G Z+, n2 G Z+ U {0}, na G Z, 
T2^(o,o,^) for ^3 G 3Z, 

where = e^*^/^ is a cube root of unity and it is easily checked that 

-712 2ni - n2 



(40) 



k 

kM 
kM^ 



27r 



\/3L 



'ni 2n2 - ni na 
' L ' V3L ' 



27r 



n2 - ni 



-ni - n2 na 



V3L 



(41) 
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FIG. 8: The last scattering surface seen from outside for a quarter- turn space E-^ with 



0.64 and = 1.28. We 



present the four pairs of matched circles along the ^-axis, which match with a twist of 7r/2, tt, 37r/2 and 27r, respectively. We 
recover the invariance by a translation of Lz for the last pair due to the fact that the cubic torus of size Lz — 1.28 is the 
four-fold cover of the quarter-turn space considered here. Only the Sachs- Wolfe contribution has been depicted here. 



One can check that 



C^''^Tn2,n2-ni -ns whcn 712 > Tlx 
C''3Tni-n2,ni,-n3 whcn m > 712 



(42) 



It follows that the analog of Eq. is given by 
1. 



n2 ,n2 — ni , — ns 



n\—n'2 ,ni, — n3 



whenn2 > ni 
when Til > 712 



(43) 



4- Sixth turn space 



The sixth turn space is like the third turn space, but with a one-sixth turn corkscrew motion 
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(44) 
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FIG. 9: The last scattering surface seen from outside for a third-turn space E4 with Lx — Ly — 0.64 and Lz — 1.92 in units of 
the last scattering surface. The first column presents the three pairs of circles alonq the z axis which match after a rotation of 
27r/3, 47r/3 and 2ti respectively. We recover the invariance by a translation of Lz due to the fact that the cubic torus of size 
Lz — 1.92 is the three- fold cover of the third-turn space considered here. The second column present the three pairs of circles 
related by translations in the x^-plane. Only the Sachs- Wolfe contribution has been depicted here. 
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The same reasoning as before shows the eigenbasis, ^^^^,^1^ k^-> 



^ [Tk + C^TkM + C'^^TkM^ + C'"^TkM3 + C'"^TkM4 + C'^^TkM^] for m G Z+, n2 G Z+ U {0}, n2 < ni,n3 G 
T2^(o,o,5^) ^3 ^ 

(45) 

where ^ = e^*^/^ is a sixth root of unity and it is easily checked that 



k 

kM 



27r 
27r 



-712 2ni - 712 ^3 



ni - n2 ni + n2 ns 



kM^ 

kM^ = 27r ( ^, 
kM^ = 27r 
kM^ 



ni 2n2 - ni na 



27r 



n2 n2 - 2ni ns 
n2 - ni -ni - n2 ns 
— ni ni — 2n2 



One can check that 



(46) 



(47) 



SO that when 723 7^ an analog relation exists with the corresponding complex random variables, and when na = 0, 
the random variable eni,n2,o is real. 



5. Hantzsche-Wendt space 



The fundamental polyhedron of the Hantzsche-Wendt space is a rhombic dodecahedron circumscribed about a 
rectangular box of size Ly/2,Lz/2). The holonomy group is generated by the three half-turn corkscrew motions 
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(48) 



The composition of these three generators is the identity, so any two suffice to generate the group. Each element 
of the Hantzsche-Wendt group has a rotational component M G {diag(l, 1, 1), diag(l, — 1, — 1), diag(— 1, 1, — 1), 
diag(— 1, —1, 1)}. The pure translations (elements with M = diag(l, 1, 1)) form a subgroup of index 4; the 
corresponding four- fold cover of the Hantzsche-Wendt space is a rectangular 3-torus of size (L^^, L^y, L^), whose 
holonomy is generated by the squares of the above corkscrew motions. Thus we may begin with the eigenspace 
for a rectangular 3-torus and ask what subspace remains fixed by the three Hantzsche-Wendt generators 

The Invariance Lemma shows that all three generators preserve the (typically four-dimensional) subspace 



^/e^} as a set, and fix the linear combination 



T/c,,/c^,/c, + (-l)""" ""^^/cx -/c^ -/c, + (-l)""^ """^-k^.ky-k, + (-1)'''" 



(49) 
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FIG. 10: The last scattering surface seen from outside for a sixth-turn space with Lx — Ly — 0.64 and Lz — 1.92 in units of 
the last scattering surface. We presents the six pairs of circles along the z axis which match after a rotation of 7r/3, 27r/3 . . . 27r 
respectively. We recover the invariance by a translation of Lz due to the fact that the cubic torus of size Lz — 1.92 is the 
six-fold cover of the sixth-turn space considered here. Only the Sachs- Wolfe contribution has been depicted here. 
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Visualizing the four wave vectors in the subscripts of as alternate corners of the cube {±kx, ^ky, ±^2) = 



27r(^^, ^^), one sees that subspace will be degenerate if and only if at least two of the indices {nx^ny^riz} 

are zero. In the degenerate case a two-dimensional subspace like {T/e^^co, "lf"-fea3,o,o} is preserved as a set, while the 
mode T/c^^o,o + '^-/c^,o,o is preserved by all three generators if and only if Ux is even. Thus the Hantzsche-Wendt's 
eigenspace has basis, t{^^^|. , 



for {fix^riy e Z+, riz G Z) or {rix = 0, riy, G Z+) or {riy = 0, rix, G 



1 

1 

1 
71 



^27r(£j,0,0) + ^27r(-^,0,0) 
^27r(0,ij,0) +^27r(0,-i^,0) 
^27r(0,0,5^) + T27r(0,0,- J^) 



One can easily check that 



for Ux G 2Z^, 
for Uy G 2Z+, 
for G 2Z+. 



2 ^[E^e] 



(50) 



(51) 



when (ux^ Uy G Z+, G Z) or (n^^ = O^Uy^Uz G Z+) or (n^y = 0, n^^, G Z+) and that otherwise ^^ffi^ is real. It 
follows that the analog of Eq. is given by 



1. when Ux^riy G and G Z, ek is a complex random variable satisfying 

^K,ky,kz — ( — 1)^"^ ^^^ka,,ky,-k^' 



(52) 



It is thus a real random variable if /c^ = and n^^ G 2Z and a purely imaginary random variable if /c^ = and 
Tlx ^ 2Z 



2. when Ux = and Uy^Uz G ek is a random variable satisfying 

so that it is a real random variable when Uy G 2Z+ and purely imaginary otherwise, 

3. when Uy = and rix^riz G Z+, ek is a random variable satisfying 

so that it is a real random variable when G 2Z+ and purely imaginary otherwise, 

4. when n^; G G Z+ or G e/c^^co, ^o,ky,o and eo,o,/c;s are real random variables. 



(53) 



(54) 



6. Klein space 



Klein space is generated by two glide reflections 
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along with a simple translation 
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(55) 
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FIG. 11: The last scattering surface seen from outside for a Hantzsche-Wendt space with Lx — Ly — Lz — 0.64 in units of 
the last scattering surface. The geometry of the circle pairings corresponds to the geometry of the fundamental polyhedron. 
The pairs of matching circles are shown here looking from directions parallel to the x-, y-^ and z-axes, respectively. Only the 
Sachs- Wolfe contribution has been depicted here. 
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The first (resp. second) ghde reflection corresponds to the two upper (resp. lower) faces of the hexagonal prism in 
Figure Q taking one to the other so that the small dark-colored (resp. light-colored) windows match. The simple 
translation takes the front hexagonal face to the back hexagonal face so that the doors match. 

The square of either glide reflection is a horizontal translation T = (La^,0,0) (taking the unmarked left wall to 
the unmarked right wall in Figure^, while the composition of one glide reflection with the inverse of the other is 
a vertical translation T = (0,L^,0). Thus the Klein space is the two-fold quotient of a rectangular 3-torus of size 

To flnd the Klein space's eigenmodes, we begin with the modes (pi?|) of the rectangular 3-torus and ask which 
remain invariant under the glide reflections (|55|) . The Invariance Lemma shows the Klein space's eigenmodes have 
the orthonormal basis 



V2 [^27r(5j,5j,5^) 



( _A ^nx+ny ^ 



for Uy G Ux^riz G Z, 
for Ux G 2Z, Uz G Z. 



(57) 



One can easily check that 



^ k k k 



(-1) 



[El 



~ kx ) ky , k z 



(58) 



when Uy G Z+, Ux^riz G Z and that ^^^ji^ j.^ is real otherwise. It follows that the analog of Eq. (^5)) is given by 



1. when Uy G Z+, Ux^riz G Z, ek must satisfy 



(59) 



It is thus a real random variable when kx = kz = and Uy G 2Z and a purely imaginary random variable when 
when kx = kz = and Uy ^ 2Z, 



2. when n^^ G 2Z, = and G e/^^ o kz is a real random variable. 



7. Klein space with horizontal flip 



The Klein space with horizontal flip is a two- fold quotient of the plain Klein space. It includes the same two glide 
reflections ((55|) as the Klein space, but adds a square root 



(60) 



of the Klein space's L^-translation Because the Klein space with horizontal flip is a quotient of the plain Klein 
space, every eigenmode of the former is automatically an eigenmode of the latter (recall the reasoning of the flrst 
paragraph of Section Hv)) . Thus our task is to decide which of the Klein space's eigenmodes (|F7)) are preserved by the 
new generator (|60j) . The Invariance Lemma shows the orthonormal basis to be 
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1/ 




V J 



(-1) 



Lix J^y z ' ^ J^x J^y z 
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1 



V L X '' L y ^ L z 



for 
for 
for 
for 



-27r(0,;t^,^) ^ ^-L27r(0 

' Ly ' Lz » ~ 

T27,(^ o,f^) + (-1)'^"T27,(_^ o,f^) 

^ J-i X z X z 

-27r(0,0,5f) 

Following the same procedure as before, we obtain that the analog of Eq. (|25|) is given by 
1. when Ux^riy G G Z, ek must satisfy 

^kx,ky,kz ~ V ^) ^kx,ky,-kz' 



'^xi'^y ^ Z^ ^ Uz G Z^ 

Uy G Z+, Uz G 2Z, 



(61) 



Tlx e 2Z+, Uz e z, 
n, e 2Z. 



(62) 



It is thus a real random variable when Uz = and Ux -\-ny G 2Z and a purely imaginary random variable when 
Uz = and rix -\- riy ^ 2Z, 
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FIG. 12: The last scattering surface seen from outside for a Klein space space Er with Lx = Ly = Lz =0.64. We present 
four pairs of matched circles. The lower right figure corresponds to the translation along the z axis while the three other plots 
corresponds to the holonomies acting in the xy-plane. Only the Sachs- Wolfe contribution has been depicted here. 

2. when Ux = 0, n^y G Z+ and Uz G Z, eo^ky,k^ must satisfy 

eS,fe„fc. = (-ir''eo,fe„-fc., (63) 

SO that it is a real random variable when = and Uy G 2Z and a purely imaginary random variable when 
riz = and Uy ^ 2Z, 

3. when Ux G Uy = and Uz G Z, e/c^,o,/c^ must satisfy 

eL.o.fe. =(-l)"^efe,A-fc., (64) 
SO that it is a real random variable when /c^ = 0, 

4. when G 2Z it has to satisfy 

eS,o,/c, = eo,o -fe.. (65) 
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FIG. 13: The last scattering surface seen from outside for a Klein space with horizontal flip Eg with Lx — Ly — Lz — 0.64. 
We present four pairs of matched circles. The first three are identical to those of Klein space, while the last one corresponds 
to the transformation along the z-axis (the upper and lower circles are mirror images of each other). Only the Sachs- Wolfe 
contribution has been depicted here. 



8. Klein space with vertical flip 



The Klein space with vertical flip replaces ((Sn|) with 




(66) 



whose action interchanges the modes ^{k^,k ,k^) ^ (— ,k^)' Consistency with the glide reflections 
whose action interchanges T(/.^ /.^ /g^) ^ {—'^)^''~^^^^{k^,-ky,k^)i requires n^c + = riz (mod 2). Thus the orthonormal 
basis is 



1 

V2 



- 27r( 



L X ^ Li y ^ L z ^ 



for Uy G Ux^riz G Z, 
for rix^riz G 2Z. 



: Uz (mod 2), 



(67) 



Following the same procedure as before, we obtain that the analog of Eq. (|25|) is given by 



1. when Tiy G and Ux^Uz G Z with 



: Uz (mod 2), ek must satisfy 



^kx,ky,kz ~ ( ^)^^^-kx,ky,-kz 



(68) 
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so that it is a real random variable when kx = = 0^ 
2. when Ux^riz G 2Z, it has to be such that 

^L,o,/c. = e_/e^,o -/c.. (69) 



9. Klein space with half turn 
The Klein space with half turn replaces or with 

f x\ /-I 0\ f x\ ( 0\ 



y 



0-1 

V ly 



y 





V J 



(70) 



The or t ho normal eigenbasis, which differs only slightly from that of the Klein space with horizontal flip (|61|) . is 



H-ir^^2H-^,-^,i^) + (-i)-^+-^+-^T,^(_ 



^27r(0,^,^) + (-l)'^''^27r(0,-^,^) 

Ivy 1^2; -f^y -f^ z 

^27r(ij,0,^^) + (-1)''^T2^(-5J,0,^) 



1 

V2 

1 
72 

^27r(0,0,Jj) 

Following the same procedure as before, we obtain that the analog of Eq. (|25|) is given by 

1. when Tix^Tiy G and G Z, ek must satisfy 

SO that it is a real random variable when /c^ = 0, 

2. when riy G Z+, G Z and Uy = riz (mod 2), ek must satisfy 



for Tlx, fly G G Z, 

for n^y G Z+, e Uy = Uz (mod 2), (71) 

for G 2Z+, G Z, 

for Uz G 2Z. 



(72) 



(73) 



so that it is a real random variable when kz = and Uy G 2Z and a purely imaginary random variable when 
kz = and riy ^ 2Z. 



3. when Ux G 2Z+ and Uz G Z, ek must satisfy 

SO that it is a real random variable when kz = 

4. when G 2Z, it has to be such that 



(74) 



(75) 



VI. DOUBLY PERIODIC SPACES 

We will first find the eigenmodes of the chimney space, and then use them to find the eigenmodes of its quotients. 
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FIG. 14: The last scattering surface seen from outside for a Klein space with vertical flip Eq with Lx — Ly — Lz — 0.64. We 
present the six pairs of matched circles. The first two rows correspond to the same holonomies as for the Klein space, while the 
last row corresponds to the translation plus flip specific to this space. Two different views of the same holonomy are shown. 
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FIG. 15: The last scattering surface seen from outside for a Klein space with half turn with Lx — Ly 



0.64. We 



present the six pairs of matched circles. The two upper left rows correspond to the same holonomies as for the Klein space, 
while the last row corresponds to the translation plus half turn holonomy specific to this space. Two different views of the 
same holonomy are shown. 
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A. Chimney Space 



Just as a 3-torus is the quotient of Euchdean space under the action of three hnearly independent translations 
Ti, T2 and T3, a chimney space is the quotient of by only two linearly independent translations Ti and T2. 
Its fundamental domain is an infinitely tall chimney whose cross section is a parallelogram (Figure [^J. And just as 
the allowable wave vectors k for an eigenmode Tk of a 3-torus were defined by the intersection of three families of 
parallel planes (Section IV A|) . the allowable wave vectors k for an eigenmode Tk of a chimney space are defined by 
the intersection of two families of parallel planes. Thus the allowable wave vectors form a latticework of parallel lines. 

The most important special case is the rectangular chimney space generated by two orthogonal translations 



Ti 
T2 

in which case the allowed wave vectors k take the form 



(ix,0,0), 

: (0,Ly,0), 



= 27r I — , 



Hx riy 

for integer values of Ux and rty and real values of rz- The corresponding orthonormal basis 



^-2'k(^ ^ ^x.riy e Z, r^ e R 



(76) 



(77) 



(78) 



is continuous, not discrete as for the 3-torus. Nevertheless, restricting to a fixed modulus /c = |k| recovers a finite- 
dimensional basis. 

The next four spaces are quotients of the chimney space, so their eigenmodes will form subspaces of those of the 
chimney space itself. 

As in the case of the torus, the random variable ek is a complex random variable satisfying 



ek = e_k. 



(79) 



B. Quotients of the Chimney Space 

1 . Chimney space with half turn 



The chimney space with half turn (Figure [^J is generated by the rectangular chimney space's x-translation Ti 
{Lx^ 0, 0) along with 
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\ Ly/2 j 



(80) 



Even though the eigenmodes are not discrete, the Invariance Lemma applies exactly as in the compact case, giving 
the eigenbasis 



1 

72 



^27.(1^ r ) + (-l)''^^27rr-Ii^ ^ -r ) i^x ^ G R) OT {Ux = 0, G 



for riy G 2Z. 



(81) 



"27r(0,^,0) 

This case is analogous to the case of the half turn space so the analog of Eq. is given by 

1. when (ux G G R) or (rix = 0, G R'^\ ek satisfies 

^K,ky,k, = {-'^T''ek^-ky,k,' 

It is thus a real random variable when Uy = and complex otherwise, 

2. when Uy G 2Z, ek satisfies 

^0,ky,0 = ^0-ky,0' 



(82) 



(83) 
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2. Chimney space with vertical flip 
The chimney space with vertical flip (Figure is generated by the translation Ti = (1/2^,0,0) along with 







(l 












1 













(84) 



The Invariance Lemma shows the orthonormal eigenbasis to be 



for Ux G Uy G 2Z. 



(85) 



One might also consider a chimney space with a vertical flip in the x-direction as well as the ^-direction. Surprisingly, 
such a space turns out to be equivalent to a chimney space with a single flip, but with cross section a parallelogram 
rather than a rectangle. In other words, the chimney space with two flips has the same topology as the chimney space 
with one flip, even though they may differ geometrically. 

Concerning the properties of the random variable, the analog of Eq. (^5)) is given by 

1. when Ux^riy G Z and G i?^, ek satisfies 

it is thus a real random variable when kx = ky = 0, 

2. when Ux ^ Z and Uy G 2Z, ek satisfies 



(87) 



3. Chimney space with horizontal flip 
The chimney space with horizontal flip (Figure (21) is generated by the translation Ti = (La^,0,0) along with 

f x\ /-I 0\ f x\ 



y 



1 

V ly 

The Invariance Lemma shows the orthonormal eigenbasis to be 



1 

V2 



( *^ ^ 

y + Ly/2 
V / 



for Ux G Z+, fly e Z^ rz ^ 
for Uy G 2Z, G R. 



"27r(0,iJ,r.) 

Concerning the properties of the random variable, the analog of Eq. (|25|) is given by 

1. when Ux ^ Z ^ Uy G Z and G i?, ek satisfies 

^L,ky,k, = {-'^T'^^kx-ky-k,, 

it is thus a real random variable when ky = kz = 0, 

2. when Uy G 2Z and G R, ek satisfies 



(88) 



(89) 



(90) 



(91) 
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4- Chimney space with half turn and flip 
The chimney space with half turn and flip (Figure is generated by 
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(92) 



(93) 



It is a four-fold quotient of the plain chimney space, unlike the preceding examples, which were two-fold quotients. 
The Invariance Lemma gives an orthonormal basis for its eigenmodes 
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^27r(0,^,r,) + 


^27r(0,5^,-r.) 
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(94) 
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V2 


T2^(j. 0) ^ 








for 


rix 






^2^(0,^^,0) 








for 


Uy 


e 2Z. 





Following the same procedure as before, we obtain that the analog of Eq. (|25|) is given by 

1. when Tlx G Uy G Z and G i?^, ek must satisfy 

^k^,ky,k^ — (~l)^^^fc^,-/cy,fe;,, 

so that it is a real random variable when ky = 0, 

2. when Uy G 2Z and G i?^, ek must satisfy 

^o,ky,k. = eo,-/cy,/c,, 
so that it is a real random variable when ky = 0, 

3. when Ux G 2Z+ and Uy G Z, ek must satisfy 

so that it is a real random variable when ky = 0, 

4. when Uy G 2Z, it has to be such that 

^0,^,0 = ^0-ky,0- 



(95) 



(96) 



(97) 



(98) 



VII. SINGLY PERIODIC SPACES 



We will first find the eigenmodes of the slab space, and then use them to find the eigenmodes of the slab space with 
flip. 
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A. Slab Space 



Just as a 3-torus is the quotient of Euclidean space under the action of three linearly independent translations and 
a chimney space is the quotient by two translations, a slab space is the quotient of by a single translation. Its 
fundamental domain is an infinitely tall and wide slab (Figure with opposite faces identified straight across. The 
allowable wave vectors k for an eigenmode Tk of a slab space define a family of parallel planes. 

If we choose coordinates so that the translation takes the form 



(0,0, 



then the allowed wave vectors k are 



k = 27r ( r^, Ty, — 



for real values of rx and ry and integer values of n^. The corresponding orthonormal basis is 

^r^^ry,^) for Tx.ry e R, Hz e Z. 



(99) 



(100) 



(101) 



Even if we restrict to a fixed modulus k = |k|, the eigenmodes of slab space remain continuous, not discrete. 

If desired, one could construct a more general slab space by identifying opposite faces with a rotation. Such a space 
would be topologically the same as a standard slab space, but geometrically different. 

As in the case of the torus and of the chimney space, the random variable ek is a complex random variable satisfying 



Gk = e-k- 



(102) 



B. Slab Space with Flip 



A slab space with flip is generated by 
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(103) 



The Invariance Lemma provides the orthonormal basis for its eigenmodes 

for Tx G Ty e Uz ^ 
-L27r(o,r^,£j) for G i?, G 2Z. 

Concerning the properties of the random variable, the analog of Eq. (|25|) is given by 

1. when rx G ry e R and G Z, ek satisfies 

^k^,ky,k^ — {~^)^^ ^k^-ky-k^- 
li is thus a real random variable when ky = kz =0, 

2. when ry ^ R and G 2Z, ek satisfies 

^0,ky,k^ — ^Q-ky-kz' 



(104) 



(105) 



(106) 



VIII. NUMERICAL SIMULATIONS 



We can now compute the correlation matrix and simulate CMB maps for the 17 multi-connected spaces described 
in the previous sections. 

In all the simulations, we have considered a flat ACDM model with I^a = 0.7, a Hubble parameter Hq = 
100/ikms~^ Mpc~^ with h = 0.62, a baryon density c^b = ^b^^ = 0.019 and a spectral index ng = 1. With 
this the radius of the last scattering surface is i?LSS = 15.0 Gpc. 
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We present a series of CMB maps with a resolution of i = 120 for the different spaces. These maps are represented 
on a sphere portraying the last scattering surface seen from outside in the universal cover. Images of the last scattering 
surface under the action of a holonomy and its inverse are shown and their intersection gives a pair of matched circles. 
All the plots presented here contain only the Sachs- Wolfe contribution and omit both the Doppler and integrated 
Sachs- Wolfe contributions. 

For the compact spaces, the characteristics of the fundamental polyhedra are 



• Lx = Ly = 0.64 for all spaces, 




• Lz = 1.28 for the half- and quarter-turn spaces (Figs. [3 IH|) 

• Lz — 1.92 for the third- and sixth-turn spaces f Figs. 151 and [Tn|) 

• Lz = 0.64 for the Hantzsche-Wendt and Klein spaces fFigs. ITTl through n3|) . 

Turning to the C^, let us examine the effects of the topology and of the volume of the fundamental domain. 

To understand the properties of the angular power spectrum on large scales, let us develop a simple geometrical 
argument based on the properties of the eigenmodes of the Laplacian operator (see Ref. j^l^ for an analogous discussion 
and Ref. 10] for the spherical case). In the simply connected Euclidean space E^, the number A^sc of modes between 
k and k + A/c is simply given by N^dk) = 47r/c^A/c, whatever the scale. Now, due to the topology, most modes will 
disappear from the spectrum and we are left with wavenumbers of modulus 

(107) 

On very small scales (large /c), the Weyl formula 8] allows us to determine the number A^^c of modes remaining 
in the spectrum: asymptotically, A^mc(< k) ~ Vk^/67r'^ (see, e.g.. Fig. 2 of Ref. fH). It follows that the number 
of modes between k and k + Ak is now given by A^mc,oo(^) ^ F/c^A/c/27r^ ~ VA'sc/(27r)^. Thus we may set the 
overall normalization on small scales where the effect of the topology reduces to an overall rescaling. But this has 
implications concerning the large scales. 

Consider a rectangular torus with a square cross section of size = Ly and with height L^, and let the relative 
proportions of Lx and Lz vary. 

When Lx ^ Lz, the space looks like a slab space and the modes on large scales (i.e., such that 27r/Lx <C k <C 27r /Lz) 

have a modulus k ^ ^ \J'^x + '^y^ ^^^^ they approach a two-dimensional distribution. Since the number of modes 

with k < ko is given by A'(< ko) = Xln=o^2(^) = tt/cq + O{ko), where rn{p) is the number of representations of p 
by n squares, allowing zeros and distinguishing signs and order (e.g., r2(5) = 8 and r3(4) = 6), we obtain that the 
number of modes between k and k + Ak is now given by A'mc,o(^) ^ LlkAk/27r. Defining the relative weight as 

(2^)^ A^Mc,o(fe) 

"^'^ = ^^Vs;(^' ^'''^ 

we obtain that w ~ (tt/ZcL^) 1 so that the large scale modes are boosted compared with the mode distribution of 
the simply connected space exactly as if the spectral index rig were lowered by 1. In the hypothesis of a scale invariant 
spectrum ng = 1, one therefore expects that the £{£ -\- 1)C£ spectrum will behave as for the relevant scales 

When Lx <C Lz, the space looks like a chimney space and the modes on large scales (i.e., such that 27r /Lz <C 
k <C 27t/Lx) have a modulus k ~ 27rnz/Lz so that they approach a one-dimensional distribution. It follows that the 
number of modes between k and k-\- Ak is now given by A'mc,o(^) ^ LzAk/n, so that w{k) ~ (27r//c^L^) 1. Again, 
this will imply a relative boost of the spectrum on large scales as if the spectral index were lowered by 2. 

When Lx ^ Lz, as long as we are above the mode cut-off, one has a three-dimensional distribution of modes so that 
the relative weight of large scale modes is ~ 1, as in a simply connected space. The signature of the topology in 
the Ci exists at sufficiently large scales in the form of small spikes around the expected value in a simply connected 
space due to the discrete nature of the /c-spectrum. 

These results are summarized in Fig. ^| 



IX. LOCATION OF THE OBSERVER 



The 3-torus, chimney space, and slab space are exceptional because they are globally homogeneous. A globally 
homogeneous space looks the same to all observers within it; that is, a global isometry will take any point to any other 
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FIG. 16: CMB anisotropies in a rectangular torus with Lx = Ly and Lz < Lx (left) of Lz > Lx (right). The volume of 
the fundamental domain is always the same (0.262i?Lss)- The spectra start from a cubic torus and Lz varies by a factor 2 
from each spectrum to the next (so that the ratio Lz/Lx reaches 64^^ for the most anisotropic configurations showed here. 
The spectrum is boosted and behaves as of depending on whether the fundamental domain is flattened or elongated 
along the z direction. Note that the scale at which the cutoff in the spectrum occurs increases with max(Lx, Lz) and that the 
spikes in the spectrum are less present for large Lz. This is because the mode spacing decays as min(L~^, L~^) and hence the 
discreteness of the spectrum is less obvious. This is of course all the more true when there are two large directions where the 
mode density is higher (left panel). 




FIG. 17: The repeating images of the last scattering circle in a two-dimensional Klein bottle align in rows if the observer 
happens to sit on an axis of glide symmetry (left), but form a different pattern if the observer sits elsewhere (right). 

point. The remaining multi-connected flat spaces, by contrast, are not globally homogeneous and may look different 
to different observers. For ease of illustration, consider the two-dimensional Klein bottle: the self-intersections of the 
"last scattering circle" are different for an observer sitting on an axis of glide symmetry (Figure El left) than for 
an observer sitting elsewhere (Figure [T71 right). Analogously in three dimensions, the lattice of images of the last 
scattering surface may differ tremendously for observers sitting at different locations within the same space. The 
power spectrum, the statistical anisotropies, and the matching circles may all differ. 

Moving the observer to a new basepoint would needlessly complicate existing computer software for simulating 
CMB maps. It is much easier to move the whole universe, leaving the observer fixed! In technical terms, we want 
to replace an eigenmode Tk(x) with the translated mode Tk(x + Xobs), where Xobs is the desired location for the 
observer. The translated mode is quite easy to compute: 

Tk(x) ^ Tk(x + Xobs) 

— gik-(x+Xobs) 

= e^^-"°^^ Tk(x). (109) 

For a simple mode Tk(x), the translation produces a phase shift (by a factor of o^'^-^ohs^ Q^^id nothing more. The full 
effect is seen when one considers linear combinations of simple modes: 

aiTk,(x)+a2Tk,(x) ^ aie^'^-''°'==Tk,(x) + a2e''^-''°'==Tk,(x). (110) 
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FIG. 18: CMB anisotropies in a half-turn space with Lx — Ly — 0.64, Lz — 1.28 for various positions of the observer. The 
observer starts from the center of the fundamental domain and moves along the x axis. As the position and size of some of 
the matching circles vary, the isotropic part of the angular spectrum also varies. The global structure of the spectrum remains 
unchanged, but the local "spikes" in the spectrum which originate from the discrete nature of the /c-spectrum are more or less 
smoothed depending on the position of the observer. For a better visibility, each spectrum has been offset by a factor of 2 
relative to the preceding one (vertical units are arbitrary). 

Each term undergoes a different phase shift, so the final sum may be qualitatively different from the original. 

Note that the phase shift (|110|) induced by the change of the position of the observer does not influence the properties 
of the statistical variable ek, but does influence the way a given mode contributes to a given angular scale. This is 
depicted in Figure El where the angular power spectrum is shown in a half-turn space for various positions of the 
observer. Corresponding examples of maps are shown in Fig. [T^ 

X. CONCLUSIONS 

This article has presented the tools required to compute CMB maps for all multi-connected fiat spaces. We gave 
for each space 

• the polyhedron and holonomy group, 

• the eigenmodes of the Laplacian. 

We then presented simulated maps for all of the nine compact non homogeneous spaces. On the basis of the angular 
power spectra we compared the effect of different topologies and different configurations for a given topology. We also 
implemented the effect of an arbitrary position of the observer which yields significant effects for non-homogeneous 
spaces. We investigated this effect both on simulated maps and angular power spectra. In particular, it shows 
that generically matched circles are not back-to-back and that their relative position depends on the position of the 
observer. 

All these tools and simulations will be of great help for extending the conclusions reached on the torus and to 
investigate their genericity as well as for providing test maps for any method wishing to detect (an interpret) the 
breakdown of global isotropy. 
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FIG. 19: The last scattering surface seen from outside for a half-turn space E2 with Lx = Ly = 0.64, Lz = 1.28 as in Fig. |3 
The observer starts from position (0, 0, 0), and slowly moves in the x direction. Due to the non- homogeneity of the space, the 
CMB maps look different. 
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